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There have been several attempts to introduce approximations into the exact foe n of 
Lilley’s equation in order to express the source term as the sod of a quadrupole whose 
strength is quadratic in the fluctuating velocities and a dipole whose suengrh s 
proportional to the temperature fluctuations. The purpose of this note is to show that 
it is possible to choose the dependent (i.e., the pressure) variable so that this type of 
result can be derived directly from the Euler equations without introducing aey 
additional approximations. 

1. Introduction 


The subject of aeroacoustics was first put on a rational basis by Lighthfll (1952, 
1954) when he rearranged the Navier-Stokes (Euler) equations into the form of a 
linear wave equation for a medium at rest with a quadrupole type source teen. The 
crucial step in Lighthill’s so called acoustic analogy approach amounts to asomkg 
that the source term is in some sense known or that it can at least be modeled in sense 
approximate fashion. While this approach was remarkably successful in 
gross features of the sound radiation from turbulent air jets, the commerci al^ 
industry ultimately found that they needed a much more sensitive tool tint was 
capable of predicting how even relatively small changff in the flow would dko ^ 
radiated sound. This motivated generations of researchers to seek improvements in tke 
Ligh thill approach. Early efforts were focused on accounting for mean now- 
interaction effects and there were a number of attempts to aoccmplish this by applywg 
ad hoc corrections to the original Lighthill predictions. A more satisfying approach 
was the one adopted by Phillips (1960), Lilley (1974) and others, which amounred » 
deriving inhomogenous moving media wave equations for the sound genentku 

The dominant part of the Lighthill source tenn is quadratic in the total flow 
velocity, which can be decomposed into a mean plus a fluctuating component. Tie 
source function therefore contains terms that are both linear and quadratic in me 
fluctuating velocity components. Lilley (1974) argued that the linear terms, wirichare 
typically such larger than the quadratic quantities, do not actually radiate any sowed 
and should, therefore, not be included in the source function, since they would tend x> 
dominate over the much smaller quadratic terms which are the true sources 
Including them would cause the sound source to be co nt a min a te d by the sm all ho t 
inevitable errors resulting from the actual computation of these terms and, wood 
thereby lead to inaccurate predictions of the sound field. 

Unfortunately, the equation derived by Lilley has a complicated source ®ni 
(Colonius, Lele, and Moin, 1997) which is not of the physically expected from L^fte J 
sum of a quadrupole whose strength is quadratic in the fluctuating velocities and a 
dipole whose strength is proportional to the temperature fluctuations. There have been 
a number of attempts to obtain such a source term by introducing wens 
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approximations into Lilley’s equation. The purpose of this note is to show that it is 
possible to choose the dependent (i.e., the pressure) variable so that this type of result 
can be derived exactly without introducing any approximations. 

2. The Lilley Equation and Related Background Information 

Lilley (1974) showed that for an ideal gas the Navier-Stokes equations 
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denotes the entropy, c p denotes the specific heat at constant pressure K = cjc r denotes 
the specific heat ratio, t denotes the time, rs are Cartesian constaitfs, p 

denotes the pressure, p the density, v = {v,, v 2 , v 3 } the fluid velocity, e„, the viscous 
stress tensor, q the heat flux vector and 
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is the connective derivative, can be rearranged into the third order wave equation (see, 
for example, Goldstein 1976, p. 253) 
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where 


n 3-ln- 2 - (6) 

K Po 

c 2 = kRT = K p/p (7) 

is the squared sound speed, R is the gas constant, T is the temperature. ¥ represents 
the effects of entropy fluctuations and fluid viscosity, which are generally considered 
to be unimportant and are therefore neglected in the following discussion. 

Linearizing the velocity and thermodynamic variables about the unidirectional 
transversely sheared mean flow 
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V, o b n U(x 2 ,x 3 ), p = Po = constant, T = r o (r 2 jr 3 ) 


( 8 ) 


and moving the nonlinear terms to the right hand side leads to the inhomogeneous 
Pridmore-Brown (1957) equation 
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is the Pridmore-Brown operator, 
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is the convective derivative based on the mean flow velocity and 


c = kRT 0 . 
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Notice that 
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when p' == p-p o « p^ i.e., when the pressure fluctuationare sroalL 

The detailed expression for T is given in Colonius, et aL, (1997). This result is 
still exact but the source term is now very complicated and even more importantly 
does not exhibit the quadrupole form originally proposed by Lighthill (1952). 
Lighthill emphasized the importance of properly exhibiting the correct muWpole order 
of the source term before introducing specific modeling assumptions for this quantity 
and Colonius et al. (1997) showed the extreme sensitivity of the predicted sound fidd 
to the detailed assumptions about the form of the source. 

Goldstein (1984) carried out a systematic second order asymptotic expansion and 
introduced a new dependent variable to show that 
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to within second order accuracy. £ 0 is defined in equation (10X the new depokfent 
variable Jt' is defined by 
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is the fluctuating sound speed. 

The source term in this equation is identical to the one that would be produced by 
an externally applied force f- {fi.fi’Jy} 311(1 is therefore properly interpreted as a 

dipole. The first term in f represents the source that would be produced by the 
fluctuating shear stress v,'v/ and can therefore be interpreted as a quadrupole. The 
remaining term is a dipole source produced by the temperature fluctuations 
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The quadrupole source scales like v’/t, where l is a characteristics length of the 
turbulence and the second term scales like T7T o (v'V-l), where TYI. is of die order of 
the turbulence Mach number squared for cold air jets, and should therefore be 
negligible compared to the first when the flow is subsonic (Morfey, Szewczyk, and 
Tester, 1978) 

Colonius et al., (1997) showed that they could accurately reproduce' the y 
numerically predicted sound field radiated from a low Mach number shear layer by 
substituting the numerically computed values for 
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and U into (14) and numerically solving the resulting linear equation for jc'. However, 
the Goldstein expansion, on which this result is based, is, at best, only locally valid, 
since nonlinear effects eventually dominate the near field disturbances and cause the 
expansion to breakdown. And since the acoustic field depends on the global solution 
to the problem, this approach does not lead to a rigorous derivation of the basic 
acoustic analogy equation. 

3, The Exact Equation 

The purpose of this note is to show that it is possible to obtain an exact 
rearrangement of the Navier-Stokes (Euler) equations that leads to a third order 
convective wave equation with a simple source term that still consists of a velocity 
quadrupole plus a fluctuating temperature dipole by introducing an appropriate 
dependent variable to represent the pressure fluctuations. 

To this end, we neglect viscous and heat conduction effects (The final result can 
easily be modified to include these effects by adding an addition term to the source 
function.) and substitute (4) into (3) to obtain 
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Then multiplying equations (1) and (2) by p'*/p, differentiating by parts and using (7) 
and (21) shows that 
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which upon introducing (1 1) and (17) can be written as 
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and 
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where we have used equation (24) to simplify equation (25). Upon introducing the 
new dependent variables 
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and using (12) and (18), these become the inhomogeneous linearized Euler equations 
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where the externally applied force / is now given by 



(30) 

These are identical in form to the linearized equations discussed in Chapter 1 
of Goldstein (1976), where it is shown (by taking the convective derivative of die first 
equations and the divergence of the second, subtracting the results and then using the 
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second equations with i = 1 to eliminate the velocity fluctuation <n the left hand side) 
that they can be rearranged into the inhomogeneous Pndmore-Jtajwn equation 
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which is identical to (14) but with the pressure fluctuation % now given by (27) and 
the externally applied force f now given by (30) rather than by (16). Notice that 


n— t^p'lpo when p' s p~pb~*0. 
Finally, it is worth noting that (3 1) can be written as 
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where L o is the same as (10) but with c 2 replaced by c 2 = c 1 +c~ and 

/ 2-^-(l + Jt)v;v'- (34) 

dXj 

is now a pure quadrupole source, which suggests that the dipole temperature source in 
(3 1) could also be interpreted as a weak nonlinear propagation. 

Aside from the definition of the pressure fluctuation* the only difference 
between equations (14) and (31) is the appearance of the pressure fluctuation factor 
(1 4- 7 c) in the quadrupole strength (1 + 7i) v'V.. Since K should be of the order of the 
turbulent intensity squared (which is typically small compared to unity) and since 
substituting exact values of U and f into (3 1) should yield exacxfy die same result as 
the direct numerical solution for the sound field, this explains why Coloneus et al. 
(1997) were able to obtain such good agreement using the two (Efferent approaches. 

Lighthill indicated that the basis of his acoustic analogy is the demonstration 
that there is an exact analogy between the density fluctuations in any real flow and 
those produced by a quadrupole source in an ideal stationary acoustic medium. The 
present result shows that, aside from viscous and heat condurtkra effects, there is an 
exact analogy between the (p/p f K fluctuations in any real flow end the corresponding 
linear fluctuations in this quantity produced by a quadrupete phis a temperature 
dipole source in an arbitrary ideal transversely sheared mean flew. 
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